Nondiffracting beams for vortex tomography 
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We propose a reconstruction of vortex beams based on the implementation of 
quadratic transformations in the orbital angular momentum. The information is en- 
coded in a superposition of Bessel-like nondiffracting beams. The measurement of 
the angular probability distribution at different positions allows for the reconstruction 
of the Wigner function. 
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In recent years, considerable attention has been paid to 
optical vortices; i.e., beams of light whose phase varies in 
a corkscrew-like manner along the direction of propagation. 
The vortex itself is a topological point defect in the wavefront 
that manifests as a null amplitude because the phase there is 
indeterminate 111 [2] • In addition, these beams carry orbital 
angular momentum (0AM) and have a topological charge £ 
that is specified by the total phase change 27r^ along a contour 
around the vortex center. 

Apart from their fundamental importance, vortex beams 
have found numerous applications fSy. Besides, entangled 
photons prepared in a superposition of states bearing a well- 
defined 0AM provide access to multidimensional entangle- 
ment. This is of considerable importance in quantum infor- 
mation and cryptography |0, S • 

Determining the 0AM state of a vortex requires knowledge 
of the phase distribution around the singularity. Since direct 
measurement of this phase in the visible is not possible, one 
needs to rely on interferometric techniques that work reason- 
ably well both at the single-photon level 1 6] and for classical 
beams iHli]. 

The goal of this Letter is to bring forward a feasible direct 
reconstruction of an optical vortex. We recall that efficient 
methods of state reconstruction are of the greatest relevance. 
Since the first theoretical proposals, this discipline has wit- 
nessed a significant growth and laboratory demonstrations of 
state tomography span a broad range of systems ||9L[TqIi. 

The canonically conjugate variable to 0AM is the angle, 
so we shall need some basic concepts about this quantity. 
Through all this Letter we use a quantum notation, although 
the translation to a purely classical language is direct. Let us 
denote by L the 0AM along the direction z of propagation. 
For our purposes, the simplest choice is to use the complex 
exponential of the angle E = e~^^, which satisfies the com- 
mutation relation [E.,L] = E um . The action of E on the 
0AM eigenstates i^ E\t) = |^ — 1), and it possesses then a 
simple implementation by means of phase mask removing a 
charge +1 from a vortex state III2I1 . Since the integer £ runs 
from —00 to +00, ^ is a unitary operator whose eigenvectors 
10) = ^== Xl^^^e"*^'^!^) are states with a well-defined an- 
gle. A point particle is necessarily located at a single value 
of the periodic angular coordinate. However, for extended ob- 



jects, such as fields, the eigenvalue (j) might be not directly 
related to an angular position in the transverse plane, as we 
shall see. 

This canonical pair describes the physics of a planar rotor. 
The phase space is then the discrete cylinder Si x Z, where 
Si is the unit circle (associated to the angle) and the integers 
Z translate the discreteness of the 0AM III3I1 . To properly 
represent a state described by the density matrix g, we need 
to introduce a Wigner function. Following our approach in 
Ref. III4I1 , we define the Wigner function for this pair as 



mA0) = 



1 

2^ 



d^'e^'^' {^-^'/2\q\^^^'/2), (1) 



whose analogy with the standard one for position and momen- 
tum is more than evident. It is advantageous to rewrite Eq. ^ 
as 

^(^. ^) = ^Y. f '^^' e'^''*-'*'^ Q{t\ 4>') , (2) 
where g{£\ cj)') are the Fourier coefficients oiW{£^(j)), given 



by^(A(/>) 



' Y,^, e'^ ^ {£'\q\£' - £). For £ = this is 



simply the Fourier transform of the 0AM spectrum {£\ q\£). 

As one step further, assume we are able to measure L^ 
transformations on the input state followed by angular pro- 
jections, that is. 



L0{<P,0 



I exp(iCLV2) Q exp(-iCLV2)|0) • 0) 



Together with the OAM spectrum, these tomograms w((/), C,) 
provide complete information. Indeed, one can check that 



q{£,^) 



1 

2^ 



d(l)'e-'^'^' uj{(l)', (})/£), 



(4) 



so the measurement of a;(^, £) allows for the determination of 
q{£ •,(})) and hence the full reconstruction of the Wigner func- 
tion via Eq. (O. 

The feasibility of this scheme relies on the implementa- 
tion of the L^ operation, which is precisely a free rotor. Our 
starting point is a monochromatic light beam whose com- 
plex scalar profile can be written in cylindrical coordinates 



as ^(r, (/?,z) = ilj{r,(p^z) ex.p{—ikz), where the envelope 
ilj{r^(p^z) is a slowly varying function of z. In further con- 
sideration, we restrict ourselves to the paraxial regime and ne- 
glect polarization effects. 

In the angle representation, L acts as —id^p and the 0AM 
eigenstates (which constitute an appropriate basis) are conve- 
niently expressed (in the plane 2: = 0) by 



V^£(r, (p) = (r, (f\£) = Ai{r) exp{ii(p) , 



(5) 



where A^(r) is an arbitrary radial profile. Please, notice care- 
fully that the polar angle (p in the transverse plane should not 
be confused with the rotor eigenstate 0. Next, we look at 
beams for which free propagation induces an L^ transforma- 
tion, as for them the tomography will be enormously facili- 
tated. We thus require 



tl;i{r,(p,z) 



iar 



i^i{r,(p). 



(6) 



In this case the free propagation of the optical field corre- 
sponds to free evolution of the rotor. Since the free evolution 
of paraxial optical beams is determined by the wave equation 

^±i^i — 2ikdzil^i = 0, where A_\_ is the transverse Laplacian, 
we obtain 



,2fMl^^dMr) 



dr'^ 



dr 



+ (£V^ -2aM^)Ai{r) =0, (7) 



whose regular solution is given in terms of the Bessel function 
Ai{r) oc Ji{V2ak£r). Consequently, the beams we look for 
turn out to have an envelope 



'^^(r, (f, z) = Ni Ji{V2akir) e 



iar 



(8) 



As it happens with plane waves, these beams require infinite 
energy, so their practical generation relies in a spatial apodiza- 
tion of some kind. We take this into account in the normaliza- 
tion factor Nf = [27r J^ dr r Ji{V2ak ir)]-^ , R being the 
radius of this apodizing element. 

Taking a ~ 1 guarantees that we are safely in the parax- 
ial region \a\ <C 2k/ f'. In that case, the Bessel beams in- 
volved have central maxima of the size of a few millimeters 
and quadratic phases aiP ~ 10 rad/m are induced in the free 
propagation. This confirms the feasibility of our proposal. 

We can also understand the form of the vortices ^ in 
an alternative way. The paraxial free- space transfer func- 
tion is H{p^ 6) (X ex.-p{—i7rXzp'^), where (p, 0) are polar co- 
ordinates in the spatial-frequency space and A is the wave- 
length. Let us assume that all the nonvanishing frequency 
components of our field lie on a circle of radius po, that 
is, V^(p,6>) oc ^{0)5{p - po), where V^(p,6>) = J^[%lj(r,Lp)\ 
denotes the Fourier transform. Then, free propagation con- 
tributes just by an unessential overall phase factor, V^'(p, 0) = 
H{p^O)^l){p^O) = exp(— zttA^Pq) '0(^), so the transverse 
intensity distribution becomes independent of the distance 
z (and hence represents a nondiffracting beam). The par- 
ti cular case of the family (O is obtained by taking po = 
V'2^K|/(27r). 

A general input state can be encoded into an optical beam 
by creating a superposition of these basis waves ?/^(r, (p) = 



"^^i Q V^^(^7 ^)' which can be easily accomplished by an am- 
plitude spatial light modulator displaying a hologram com- 
puted as an interference pattern of the required field ?/^(r, (p) 
and a reference plane wave. 

After propagating a distance z the beam will evolve into 
a new superposition ?/;(r, (p^ z) = ^^ CiC^^^ ^ V^^(^, ^)- But 
according to Eq. (O, this is just the required L^ action with 
C = —az. This is the main result of this paper: by chos- 
ing a particular form ([8]) of the 0AM eigenstates, the highly 
nontrivial L^ transformation is realized by a simple free beam 
propagation. This also provides an experimentally feasible 
optical realization of a free rotor. 

After some simple calculations, the associated tomograms 
for the superposition field turns out to be 



•^^-h 



^Qe^^^e-<^'/2 



(9) 



By scanning the values of C, the beam can thus be recon- 
structed according to Eq. (|4]). This completes our proposed 
scheme. Naturally, Q are independent of the rapid oscillat- 
ing factor exp(— i/cz). 

As a last experimental detail, we need to assess the mea- 
surement of the angular spectrum in Eq. ([3]). For the field 
?/^(r, (/?) this reads as (taking for simplicity the plane z = 0) 



uj{4>,^) = P{4>) 



I 



dipdrriljl{r,(p)ilj{r,ip) 



(10) 



where ^^^(r, (p) = (r, (^|0), and |^) is a field with well-defined 
angle. Note in passing that the angular spectrum P((/>) cannot 
be obtained, in general, by measuring the angular profile of 
the output beam, e.g., using a wedge-shaped aperture. In fact, 
as a consequence of the nontrivial radial profile of the 0AM 
eigenstates, the states of ^^^(r, (p) are given by 



M^.^) = ^e'^ 



(11) 



Finally, using the cross-correlation theorem of Fourier anal- 
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ysis, one gets P(0) = 



T- 



. This re- 



r=0 



suit provides a simple interpretation of the tomographic setup. 
The angular states (fTTI) are defined in such a way that by plac- 
ing a transparency ^^^(p, —6) in the Fourier plane, the un- 
wrapping of azimuthal phase modulations occurs simultane- 
ously for all the vortices in the superposition ?/^(r, (p). As a 
result, the on-axis interference observed at the output plane 
of the optical processor yields phase information about the 
coefficients q. Different angular components can be readily 
accessed by rotating the filter transparency with respect to the 
optical axis. 

The proposed experiment thus consists of the following 
four steps: (i) The input state is encoded into a superposi- 
tion of Bessel-like beams V^i(r, (p). (ii) This superposition is 
propagated a distance z. (iii) The angular distribution uj{(j), z) 
is determined by measuring the intensity in the middle of 
the output plane of the optical processor as a function of the 




Fig. 1. (Color online) Calculated (left) and experimentally 
generated (right) transverse intensity distribution for a nor- 
malized superposition of the states V^i(r, (/)) and V^2(^, 0)- 



Fourier filter orientation cj). The last two steps are repeated for 
a set of distances z. (iv) The Wigner function of the input state 
is reconstructed with the help of Eqs. (O and (|4]). 

As an illustration, we apply our reconstruction procedure to 
the state '^(r, (p) = [ipi (r, (p) + '^2(^5 (p)]/V^. The calculated 
transverse intensity of this state is shown in the left panel of 
Fig.[T] while a typical CCD scan of the state generated in the 
laboratory by an amplitude spatial light modulator (CRL Opto 
1024x768 pixels) is shown on the right for comparison. As 
an illustration, we apply our reconstruction procedure to the 
state ^p{r^ip) = [?/^i(r, (^) + ^|J2{r^(f)]/V^. The calculated 
transverse intensity of this state is shown in the left panel of 
Fig.[T] while a typical CCD scan of the state generated in the 
laboratory by an amplitude spatial light modulator (CRL Opto 
1024x768 pixels) is shown on the right for comparison. 

The tomograms are calculated from Eq. (|9]), whence we 
obtain the Fourier coefficients as ^(^,0) = e*^'^^^((^^,i + 
(5^ _i)/(47r). This leads to the Wigner function 



47r 



(-1) 



^+1 



7r2(3-2£) 



COS( 



(12) 



We can see three contributions: two flat slices, coming from 
the fields i/ji and ^^2, located at £ = 1 and £ = 2, respectively. 





Fig. 2. (Color online) Reconstructed Wigner function ([T2I) 
corresponding to the same state as in Fig. 1. Left panel shows 
this function on the unit cylinder (whose vertical axis corre- 
sponds to the 0AM), while right panel show the unwrapped 
plot in Cartesian axis. 



and an interference term that contains contributions from all 
the values of £, although damped as l/£. These features are 
illustrated in Fig. [2l which shows the reconstructed Wigner 
function on the unit cylinder and its wrapped counterpart in 
Cartesian axes. 

In practice, one cannot scan infinite values of ( (i.e., of z), 
and only some discrete version of the Eq. (|4| is available in the 
laboratory. We have checked that with 50 equispaced values, 
the reconstruction is still fairly stable. 

In summary, we have carried out a full program for a com- 
plete reconstruction of vortex beams. In our scheme, the in- 
formation is encoded in a superposition of Bessel-like non- 
diffracting beams for which free propagation is equivalent to 
the L^ operation. At the same time, our formulation provides 
a theoretical framework on the action of angle and angular 
momentum on these beams. 
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